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ABSTRACT 


We present the discovery of 25 new repeating fast radio burst (FRB) sources found among 
CHIME/FRB events detected between 2019 September 30 and 2021 May 1. The sources were found 
using a new clustering algorithm that looks for multiple events co-located on the sky having similar dis- 
persion measures (DMs). The new repeaters have DMs ranging from ~220 pc cm”? to ~1700 pc cm ?, 
and include sources having exhibited as few as two bursts to as many as twelve. We report a statistically 
significant difference in both the DM and extragalactic DM (eDM) distributions between repeating and 
apparently nonrepeating sources, with repeaters having lower mean DM and eDM, and we discuss the 
implications. We find no clear bimodality between the repetition rates of repeaters and upper limits on 
repetition from apparently nonrepeating sources after correcting for sensitivity and exposure effects, 
although some active repeating sources stand out as anomalous. We measure the repeater fraction 
and find that it tends to an equilibrium of 2.655796 over our exposure thus far. We also report on 14 
more sources which are promising repeating FRB candidates and which merit follow-up observations 
for confirmation. 


Keywords: Radio transient sources (2008); High energy astrophysics (739) 


1. INTRODUCTION 


Fast radio bursts (FRBs) are microsecond—millisecond flashes of radio waves that are detectable over extragalactic 
distances (see Cordes & Chatterjee 2019; Petroff et al. 2019, for reviews of the phenomenon). Some FRB sources 
repeat, which rules out cataclysmic models for burst production (Spitler et al. 2016). Magnetars are likely to be 
the origin of at least a fraction of FRBs as the two source classes have been directly linked (Bochenek et al. 2020; 
CHIME/FRB Collaboration et al. 2020a). The sources that repeat show clustered activity (Oppermann et al. 2018) 
that in some cases is periodic, with bursts being detectable only during a window of time in each cycle (CHIME/FRB 
Collaboration et al. 2020b; Rajwade et al. 2020; Cruces et al. 2021). For one source, this periodic activity has been 
shown to be chromatic, with activity at lower frequencies trailing the activity at higher frequencies (Pastor-Marazuela 
et al. 2021; Pleunis et al. 2021a; Bethapudi et al. 2022). Importantly, it is as-yet unclear whether all FRBs repeat 
(Caleb et al. 2019; James et al. 2020a,b), what fraction of repeating sources show periodic activity, and on what 
timescales repeaters are active. 

Complex time-frequency structure, with subbursts that drift downward in frequency as time progresses, seems 
ubiquitous among bursts from repeating sources (Hessels et al. 2019; CHIME/FRB Collaboration et al. 2019a). On 
average, repeater bursts are observed to be longer in duration and narrower in bandwidth than for bursts seen from 
nonrepeating sources (Pleunis et al. 2021b). It is unclear if this observed difference represents an intrinsic difference 
in the two populations, if it is caused by a propagation effect related to distinct source environments, or if it is due to 
a selection effect that leads to the observation of two extremes out of a continuum of source properties (e.g., Connor 
et al. 2020). 

Repeating sources of FRBs provide the best opportunity for detailed follow-up observations, and the most pro- 
lific of repeaters - FRB 20121102A (Spitler et al. 2016), FRB 20180916B (CHIME/FRB Collaboration et al. 2019b), 
FRB 20190520B (Niu et al. 2022), FRB 20200120E (Bhardwaj et al. 2021a), and FRB 20201124A (Lanman et al. 2022) 
— have provided many insights. All five have been localized to sub-arcsecond precision. They are associated with a 
star-forming dwarf galaxy, a spiral, another star-forming dwarf galaxy, a globular cluster, and a barred galaxy, respec- 
tively (Tendulkar et al. 2017; Marcote et al. 2020; Niu et al. 2022; Kirsten et al. 2022; Xu et al. 2022), demonstrating 
that repeaters inhabit a variety of host types. Two sources of repeating FRBs that live in extreme magneto-ionic 
environments are coincident with persistent radio sources (PRSs; Michilli et al. 2018; Niu et al. 2022) but for other 
sources, PRSs of luminosity similar to these two have been ruled out (Law et al. 2022). The complex magnetized 
environments of repeaters can result in depolarization of bursts towards lower frequencies (Feng et al. 2022), depolar- 
ization due to Faraday conversion (Kumar et al. 2022) and temporal evolution of the Faraday rotation that can be 
secular (Hilmarsson et al. 2021), a combination of stochastic and secular (Mckinven et al. 2022a) or exhibit a sign flip 
pointing to a magnetic field reversal (Anna-Thomas et al. 2022). 
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Repeater bursts have been resolved to timescales from a few ps to a few ms (FRB 20180916B and FRB 20200120E; 
Nimmo et al. 2021, 2022a) and they have been detected up to 8 GHz (FRB 20121102A; Gajjar et al. 2018) and down 
to 110 MHz (FRB 20180916B; Pastor-Marazuela et al. 2021; Pleunis et al. 2021a). Repeater bursts can be extremely 
narrowband (65 MHz; Kumar et al. 2021) and show order-of-magnitude differences in brightness from subburst to 
subburst (Chawla et al. 2020; Caleb et al. 2020) or from burst envelope to burst envelope (Kumar et al. 2019). The 
unprecedented sensitivity of the Five-hundred-meter Aperture Spherical Telescope enables detailed measurements of 
burst energy distributions (Li et al. 2021), which informs possible progenitor models and emission mechanisms through 
the determination of minimum and maximum characteristic energies.! Despite coordinated efforts, a prompt multi- 
wavelength counterpart to a repeater burst has not been detected (e.g., Hardy et al. 2017; Scholz et al. 2020; Tavani 
et al. 2020; Niino et al. 2022). 

It has not yet been settled what observed behavior is necessary for FRB production (i.e., if an extreme magneto-ionic 
environment and/or a mechanism that provides periodic activity are a requirement or a coincidence?) and it is unclear 
whether models that have been tuned to the properties of known active repeaters can be safely extrapolated to explain 
all FRBs. A larger sample of repeating sources of FRBs is key to exploring similarities and differences between the 
population of repeaters and apparent nonrepeaters. 

Optimal strategies for discovering repeating sources of FRBs are to observe FRB sky positions with more sensitive 
instruments or to revisit the same patch of sky repeatedly (Connor & Petroff 2018). Indeed, two FRBs discovered 
by the Australian Square Kilometre Array Pathfinder were shown to repeat by the more sensitive 100-m Robert C. 
Byrd Green Bank Telescope (Kumar et al. 2019) and the 64-m Murriyang (Parkes) telescope (Kumar et al. 2021), and 
FAST discovered repeat bursts from an FRB initially discovered by the Murriyang telescope (Luo et al. 2020). 

The Canadian Hydrogen Mapping Experiment (CHIME) is a radio interferometer with no moving parts that observes 
the full Northern sky (6 > —11?) daily. Its FRB project (CHIME/FRB; CHIME/FRB Collaboration et al. 2018) has 
previously reported twenty repeating sources of FRBs (CHIME/FRB Collaboration et al. 2019a,b; Fonseca et al. 2020; 
Bhardwaj et al. 2021a; Lanman et al. 2022) through a continuous survey that commenced in mid 2018. 

As the number of detected FRBs continues to grow, the chance coincidence probability of finding two unrelated FRBs 
at similar sky position and DM becomes non-negligible. This is especially true when the sky localization precision of 
an FRB is coarse and the exposure towards a certain sky area is large. In this paper we present an improved and more 
accurate framework for the calculation of chance coincidence probabilities that builds on the approach outlined by 
Fonseca et al. (2020). We present 25 new repeating FRB sources discovered by CHIME/FRB from 2019 September 30 
to 2021 May 1. In order to systematically search our events for repeating sources, we employed a clustering algorithm 
on all events in our database. We outline our source identification procedure in 82 and source and burst characterization 
in 83. We interpret these findings in 84 and conclude in 85. 

Note that while preparing this manuscript, CHIME/FRB alerted the community to our discovery of new repeater 
FRB 202209124, which falls outside of the cutoff date for the analysis presented here (Mckinven & CHIME/FRB 
Collaboration 2022). We will therefore report on it elsewhere. 


2. OBSERVATIONS 
2.1. CHIME/FRB 


CHIME is a radio interferometer located at 49? latitude near Penticton, B.C. in Canada. It consists of four semi- 
cylindrical 100-m x 20-m reflectors with 256 dual-polarization cloverleaf antennas installed on each focal line (Bandura 
et al. 2014). The instrument is sensitive from 800 MHz down to 400 MHz. The FRB detection system (CHIME/FRB 
Collaboration et al. 2018) runs on the telescope in parallel to hydrogen mapping (CHIME Collaboration et al. 2022) 
and pulsar timing and search projects (CHIME/Pulsar Collaboration et al. 2021), and tiles the primary beam of 
the telescope with four East-West x 256 North-South stationary beams that each yield 16,384 frequency channels 
and a 0.98304-ms time resolution. These data streams are cleaned of radio frequency interference (Rafiei-Ravandi & 
Smith 2022), dedispersed with a tree algorithm using trial dispersion measures (DMs) up to approximately 13,000 pc 
cm”? and then searched for peaks, all in real time. All trigger metadata are stored in a database, but only upon 
detection of a significant and likely astrophysical trigger are data around the trigger time saved to disk from a ring 


l Li et al. (2021) claim a bimodal energy distribution for FRB 20121102A, but Aggarwal (2021) argues that the bimodality disappears when 


energies are corrected for bandwidth. 


2 A case in point is the magnetar near the Galactic center, for which the observed evolution of Faraday rotation has been demonstrated to 


be unrelated and extrinsic to the burst production (Desvignes et al. 2018). 
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buffer. CHIME/FRB has the ability to store “intensity” (the beamformed Stokes-I data streams) and “baseband” 
(the complex voltages measured at each of the 1,024 antennas) data. Whereas the former has a much smaller storage 
footprint, the latter allows for repointing the array in software (Michilli et al. 2021) and gives access to the native data 
resolution (800 x 10° samples s^, in 1,024 frequency channels at 2.56 us time resolution by default) and polarimetric 
information (Mckinven et al. 2021), which allows for coherent removal of the dispersion delay. 


2.2. Source Identification 


Our systematic identification of new sources and associated bursts is broken into four stages. First, we obtain a list 
of candidate clusters by running a clustering algorithm on the coarse sky positions and DMs from detection metadata 
in the database (82.2.1). Then, for all bursts under consideration, we refine the measurements of their sky position 
and DM using the best available data (83). These measurements then form the input for a contamination rate (Rec) 
calculation, based on the probability of chance coincidence (P«.), that determines the likelihood of nearby events being 
unrelated (82.2.2). All sources with sufficiently low Rec are investigated further. The Rec calculation is also used 
to identify any outlier events in clusters (e.g., when three events were grouped together by the clustering algorithm, 
but only two events in the cluster are from the same source). Finally, we verify the self-consistency of all combined 
localizations (83.1). All sources that pass this final test are deemed real repeaters. 

Here, we have analyzed all detection metadata in the CHIME/FRB database between 2019 September 30 and 2021 
May 1. 


2.2.1. Clustering Algorithm 


We clustered events using an in-house implementation (Dong et al. 2022) of the DBSCAN algorithm (Ester et al. 
1996). DBSCAN is an unsupervised machine learning algorithm to find clusters of points in n-dimensional space. We 
wish to find clusters in the coarse right ascension (o), declination (6) and DM and for this, DBSCAN offers two major 
benefits. First, it is insensitive to the shape of the cluster as it does nearest-neighbor clustering. Bright sources in 
CHIME/FRB are detected in arcs of apparent o and 6 due to pulses/bursts being sequentially detected in the Eastern 
sidelobes, the main lobe and the Western sidelobes of the formed beams as the sources transit through the telescope's 
field-of-view.? With DBSCAN, one such arc of detections from one source is grouped in one cluster and not in multiple 
clusters. Secondly, DBSCAN can form an unbounded number of clusters and no arbitrary human choice of number of 
clusters is required. 

Based on our knowledge of CHIME/FRB features and systematics, such as the sidelobe detections of bright sources, 
we have set the nearest-neighbor distance, e, to be a function of the measurement uncertainties for event parameters. 
The uncertainty in o is not constant across our declination range (> —11?) and increases towards the North Celestial 
Pole (NCP; ô = 90°): og = 2.2?/cos(0). The other tolerances are set to be constant fiducial values 0; = 0.57, 
opm = 13 pc cm”?. These are the approximate beam size full width at half maximum and the largest DM uncertainty 
in our tree dedispersion algorithm. The op» threshold also allows for grouping of events that were detected by the 
real-time pipeline at suboptimal DMs (see also §4.1). The minimum number of points in each cluster is set to be two. 

This method successfully recovered all previously published CHIME/FRB-discovered repeating FRBs. After the 
removal of clusters from known sources”, the analysis resulted in the identification of 65 new clusters. 


2.2.2. Chance Coincidence Probability Calculation 


As the CHIME/FRB experiment continues to detect more FRBs, the probability of identifying two or more unrelated 
FRB sources with similar DM and within a typical localization region becomes non-negligible. To calculate the 
probability that each of the x bursts in a cluster are physically unrelated to one another we treat the detection of 
FRBs as a set of independent Bernoulli trials, with each FRB detected by CHIME/FRB as one trial. For each cluster, 
we are then interested in the probability of detecting an additional x — 1 coincident bursts (a — 1 “successes” ) that are 
physically unrelated (hence, independent). Coincident here means that within some volume AaAdADM the values 
for burst parameters agree. The probability P that another coincident burst is detected can be written according to 


3 The arcs of apparent a and 6 are symmetric around, and peak in S/N, at the true a and 6 of the source. 

4 These known sources can be previously discovered repeaters or detections of known pulsars in off-meridian sidelobes, where their apparent 
sky positions makes them seem extragalactic from comparing their observed DMs to the maximum expected Galactic DM contributions 
towards their apparent sky positions. 
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the binomial distribution 


wv (N-1 k (N-1)-k 

Pec düster = 5 ( k Jetasow ( a P{c,5,DM}) , (1) 
k<z—1 
where N — 2196 is the total number of FRBs detected by CHIME/FRB? by the cut-off date of this analysis and 
P{a,5,DM} is the probability of a “success” (the probability that a source detected by CHIME/FRB is detected in the 
AaAÓADM volume). 

The probability pto,5;pw is a function of the distribution of FRB sources in the Universe and the CHIME/FRB 
exposure and selection function. We have estimated it empirically? by using a normalized histogram of our entire 
FRB source sample (i.e., without repeat bursts) smoothed by a Gaussian kernel with an angular scale of 10? and a 
DM scale of 350 pc cm^?. Ideally, we want to consider only independent FRB sources for this analysis. It is, however, 
not possible to obtain a sample completely free of repeat bursts. When repeat bursts are inadvertently left in the 
sample, that artificially inflates p{a 5, pm} in the respective portions of the sky, which leads to higher chance coincidence 
probabilities, making our analysis more conservative. 

Using N — 2196, as above, as our sample size/number of trials, we calculate a contamination rate (i.e., the expected 
number of false positives in our sample) 

Ree = Pee X N, (2) 
that is shown in Figure 1 for each cluster. We select which sources we deem to be real repeating sources of FRBs by 
setting a maximum acceptable contamination rate Rece < 0.5 (“gold” sample), such that there is a 50% probability 
that one source in the sample got included by chance. This “gold” sample includes 32 previously unpublished sources. 
We also present clusters with 0.5 < Ree < 5 (“silver” sample) in Appendix A as candidate repeaters that may warrant 
follow-up observations. This “silver” sample includes an additional 18 previously unpublished potential repeaters. 
However, upon closer inspection, seven and four of those clusters, respectively, have inconsistent event localizations 
(see 83.1 and Fig. 1). We thus end up with 25 new repeater sources (*gold" sample) and 14 new repeater candidates 
(“silver” sample). 

As it is possible that a cluster contains unrelated bursts (e.g., from the tolerances of the clustering analysis being 
large; see 82.2.1), we calculate Rec for each combination of events" in a cluster and consider the minimum Ree of all 
combinations (which should be the Rec of all events in a cluster, if no event is an outlier). We inspect all combinations 
(see Fig. 14 in Appendix E) and flag those events that increase the Rec by more than an order of magnitude as 
compared to the minimum Rec of that cluster in one or more of the combinations (ten potential outliers in total in 
nine unique candidate clusters; these flagged events are marked with asterisks in Figs. 5, 15-17). However, we do 
not exclude from our sample those flagged events, as a mix of parameter uncertainties (sky positions from real-time 
detection S/Ns versus baseband data and a S/N-dependent DM uncertainty) could also lead to a “forked path" of Rec 
as a function of included events, where one event leads to a systematically lower Rece when combined with other events 
in the cluster (see Fig. 14 for examples). The outlier analysis would also reveal if two or more separate repeaters were 
grouped into one cluster, but we see no evidence for this. 


3. ANALYSIS & RESULTS 


'The properties of 25 new repeating sources of FRBs are summarized in Table 1 and individual burst properties are 
collected in a machine-readable table?, an excerpt of which is shown in Table 5 in Appendix C. A sky distribution of 
repeating sources of FRBs is presented in Figure 2. Dedispersed dynamic spectra of the bursts, best-fit localization 
regions of each source and the exposure up to 2021 May 1 are publicly available.? 

Figure 3 shows the detection timeline of all repeating sources of FRBs discovered and detected (i.e., FRB 20171019A, 
originally discovered by ASKAP; Kumar et al. 2019) by CHIME/FRB up to 2021 May 1, including previously unpub- 
lished bursts from already published sources. Of the previously unpublished bursts from the known sources, we have 
measured their burst properties if their first detection with S/N > 10 or 12 happens to be an unpublished burst for 


5 This number is subject to change as we vet FRB candidates for future catalogs. However, the Pec calculation is only marginally sensitive 


to small variations in N. 


6 With the exception of a, for which we can assume a uniform probability distribution, as our exposure covers multiple years (see also Fonseca 


et al. 2020). 
7 That is, all unique combinations (5). where n is the number of events in each cluster and 2 € k< n. 


8 https:/ /www.chime-frb.ca/repeater. catalog 
9 https: //www.canfar.net /storage/vault /list /AstroDataCitationDOI/CISTI.CANFAR/23.0004/data 
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Figure 1. The panel on the left shows the minimum contamination rate (Rec) of all combinations of bursts in 65 clusters of 
events as a function of the average declination of the cluster. Note that the vertical axis is inverted. The gray regions indicate 
0.5 < Rec < 5 (“silver” sample; lighter region) and Rec > 5 (darker region), markers are color-coded by cluster size, and clusters 
with inconsistent localizations (which are then discarded) are indicated with a black cross. The panel on the right shows a 
histogram of the Rec. In general, clusters with more bursts and/or baseband localizations and clusters at lower declination have 
lower Rec. 


the analysis in 83.3. We use the detection S/N and times-of-arrival of these bursts for the analyses in 84.4 and 84.5. 
A machine-readable table of these bursts is available in the same format and at the same location as the table with 
bursts from the new sources. 

Results on host galaxies to the new sources will be reported by Ibik et al. (in prep.), on associated persistent radio 
emission by Ibik et al. (in prep.) and on cross-correlation with galaxy catalogs by Rafiei-Ravandi et al. (in prep.). 


3.1. Localization 


For the purposes of clustering events to obtain repeater candidates, we use the positions from the real-time detection 
pipeline (“header” localization), which are based solely on the per-beam S/N values. This localization method has 
been used for past CHIME/FRB publications; most recently CHIME/FRB Collaboration et al. (2021). For an in-depth 
description of this localization methodology, see CHIME/FRB Collaboration et al. (2019b). Briefly, the method is a 
grid-search x? minimization, where observed ratios among per-beam S/N values are compared to predictions from the 
CHIME/FRB beam model””, over a grid of possible sky locations. The method can be applied to multiple bursts from 
the same source by summing the x? values for each data set. As the best-fit sky position, we construct a 90% confidence 
region. For most sources, this forms one “island”, but for a few sources (FRB 20190226B and FRB 20190113A) two 
islands are equally likely, and we present two best-fit localizations in Table 1. 

For those bursts for which baseband data were saved to disk, we run the baseband pipeline detailed by Michilli et al. 
(2021) to generate a sky-localization. This pipeline works in two main stages. First, a crude grid of tied-array beams 
is formed to determine whether the burst occurred in a sidelobe, and to obtain an initial position estimate. Next, a 
tightly-packed grid of beams is formed around the initial position to get a refined best-estimate position. In ~20% 
of cases, the baseband pipeline does not converge, due to a combination of low S/N, dedispersion artefacts and other 
reasons. Work to debug these issues is underway, and we defer those baseband localizations to future work. We do 
still include the events in our analysis, but use their positions from the real-time detection pipeline. 


10 https://github.com/chime-frb-open-data/chime-frb- beam-model 
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Figure 2. Sky distribution (Aitoff projection) of repeating sources of FRBs discovered by CHIME/FRB (orange open diamonds 
and green open circles; this work and CHIME/FRB Collaboration et al. 2019a,b; Fonseca et al. 2020; Bhardwaj et al. 2021a; 
Lanman et al. 2022), discovered by other telescopes (gray open circles; Spitler et al. 2016; Kumar et al. 2019; Luo et al. 
2020; Kumar et al. 2021; Niu et al. 2022) and detected by CHIME/FRB (purple open squares; Josephy et al. 2019; Patel & 
CHIME/FRB Collaboration 2019). 6 < —11° is outside of CHIME/FRB’s field-of-view and shaded light gray. The gray solid 
lines in the background show the plane of the Galaxy (b = 0°) as well as lines of constant Galactic longitude 0° to 360° in steps 
of 30°. 


When one or more baseband localizations are available, a combined position is calculated using the inverse-variance 
weighted mean. The uncertainties on the combined sky positions are likely overestimated, as the systematic uncer- 
tainties are accounted for in each individual error bar (Michilli et al. 2021). We defer a more precise treatment of the 
combined systematics to future work. When no baseband localization for a source is available, a final best sky position 
for each source is calculated from combining header localizations as was done by Fonseca et al. (2020).!! 

When labelling a collection of events as a candidate repeater, localization consistency is crucial to check. The 
nature of this check depends on the available data. For candidates with multiple baseband triggers, we verify the 
positions are consistent such that a weighted average yields no outliers. We also ensure any given baseband position is 
consistent with the uncertainty region derived from the header data alone, and we verify that each position is localized 
to the correct lobe of the formed beams and not a sidelobe. Baseband positions that are contained within the 90% 
confidence region are considered fully consistent, those in close proximity (within a few arcminutes) are considered 
marginal, and those several arcminutes to a few degrees away are considered outliers. For all repeater candidates, a 
successful joint-fit header localization is the baseline requirement. This test fails for seven clusters in the “gold” and 
four clusters in the “silver” sample; see §2.2.2 and Fig. 1, and those clusters are discarded. Conceptually, the joint fit is 
similar to taking the intersection of the uncertainty regions of individual bursts. Individual uncertainty regions reflect 
the instantaneous sensitivity of the beams that resulted in the detection, as well as those that did not. A merged 
localization with (near-)zero sky area is indicative of bursts that do not originate from the same sky location. We take 
all combinations of events from a cluster and check the merged localizations. If just one event led to the (near-)zero 
sky area, we discard that event. If all events in a cluster are inconsistent with each other, we discard the full cluster. 

Figure 4 shows four source localizations that represent four possible ways in which the combined header and (com- 
bined) baseband positions can converge to the true sky position of a source. In the example of FRB 20190430C (top 
left), multiple header localizations do not greatly reduce the allowed phase space, and the baseband positions land 


11 The only difference is that here, we do not attempt to rule out regions of parameter space by comparing the spectral occupancy of the 
bursts with the spectral sensitivity predicted by our beam model. 
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Figure 3. Detection timeline for 46 repeating sources of FRBs (45 discovered by CHIME/FRB and FRB 20171019A; Kumar 
et al. 2019) from the start of CHIME/FRB operations to 2021 May 1. Different sources have different colors, detections in 
the “lower” transit are indicated by a downward-facing triangle and repeater discovery dates (when the second burst above a 
completeness threshold, here S/N > 12, was detected — not the case for all sources; see §4.5) are highlighted with an extra circle 
or triangle. All 25 newly discovered repeaters are indicated with an asterisk on the right-hand ordinate axis and the names of 
some notable sources are listed. Dotted vertical lines mark the cutoff dates for CHIME/FRB Collaboration et al. (2019b) and 
Fonseca et al. (2020). 


CHIME/FRB Discovery OF 25 REPEATERS 9 


Table 1. Properties of 25 New Repeating Sources of FRBs, Ordered by Increasing Rec (Our “Gold” Sample) 


FRB Source” Ree a? Du rid b^ DM? DMk&z»o1 DMYuwie Noursts Exposure! Completeness? Burst rate” 
(J2000) (J2000) (deg) (deg) (pccm ?) (pccm ?) (pe cm~?) (hr, upper / lower) (Jy ms) (hr 3) 
20201130A — 5.3x 107^ — 64.39(1) 7.94(1) 85.4 —29.0 287.984(7 56 69 12 15411 10.4 2.3473 oe 
20191106C' 3.1 x 1077? — 199.58(1)  43.002(9) 105.7 73.2 333.40(2) 25 21 7 66-£3 1.9 (2.54*223) x 10? 
20200619A — 43x10 ! J272.6(1) 55.56(6) 84.0 27.8  439.772(4 51 45 5 110+ 2.5 (1.3*52) x 107? 
20190804E — 5.2x 107  261.34(2)  55.069(8) 82.9 34.1  363.68(1) 43 37 6 9404 1.9 (1.52* 558) x 107? 
20190915D — 3x10 !? 11.78(3) 46.86(2) 22.2 —16.0  488.69(2) 89 88 7 92 + 11.3 (2.56* 7-38) x 107! 
20200929C' 107° 17.04(2) 18.47(1) 28.5 —44.2 413.659(4 38 29 7 56 +2 5.5 (1.45* 522) x 107! 
20201221B! — 8.3x 1077 = 124.20(3) . 48.78(2) 170.5 33.8 510.42(5) 51 46 6 90-2 34 (3.787395) x 107? 
20200809E — 2.8x10-9 20.0(1) 82.89(2) 123.9 20.1  1703.48(1 74 83 4 510 +2 / 410 +26 25/641 (2.81* 185) x 107? 
20190600C — 4.8x 1079  73.32(1)  24.068(6) 177.4 —12.4 479.602(9 113 154 3 7240 24 (7.3* 557) x 107? 
20190226B  &T7x10 5 273.6201] 61.677911 90.9 28.0  632.597(9 51 45 3 85 +52 14 (3.65*$29) x 107? 
270.2502;  61.03*01; 911 29.6 48 43 
20200223B' — 41x10 ? . 8.265(8) 28.831(7) 118.1 .-—33.9 202.268(7 46 37 5 524+ 7.6 (1.447558) x 107! 
20200202A 0.00028 25.93(2) 44.290(7) 132.7 -17.6 722.37(7) 83 84 3 6245 6.1 (6.535 152) x 107? 
20190430C 0.00047  277.210(9)  24.770(9) 53.1 15.7  400.57(1) 100 84 2 29 + 11.5 (2.36* 253) x 107! 
20200420A 0.0011 7.57*004 — 37.20*00& 118.3 —25.4 671.5(4) 58 50 2 55 + 29 0.7 (9.587363) x 1074 
20200127B 0.017 119.2(1)  86.609(8) 126.6 280 3513(4) 54 51 2 154+ 13 / 951 +18 3.7 / 5.3 1.05 x 107? 
20201114A 0.018 221.59(8)  71.79(3) 11.2 42.6 322.23(1) 38 31 2 135 +27 / 83417 6.8 / 66.3 (2:837 3:08) x 19-2 
20190110C 0.021 249.33(1) | 41.445(9) 65.6 42.1  221.92(1) 37 30 3 4511 1.2 (8127) x 107? 
20200118D 0.038 106.91(1) 42.837(9) 174.4 20.9  625.57(3) 77 92 2 723 1.2 (3.3723) x 10? 
20181226F 0.049 95.87*023  83.80*017 129.8 26.2  241.887(3) 58 56 3 394 + 192 / 376 + 215 4.8 / 10.6 (2.37* $86) x 107? 
20191013D 0.051 40.42(2) 13.63(3) 159.7 —41.3  523.57(3) 43 36 2 67+ 8.9 (7.097 533°) x 107? 
20191114A 0.066 273.57(3) 19.80(4) 47.0 16.8 552.47(6) 98 80 3 6145 21 (1.36* 2 5) x 107? 
20190113A? 0.13 108.26*020  —2.98*021 218.1 3.4  428.90(1) 178 251 2 40 + 25 5.0 (4.97* 2,39) x 107? 
106.50*9:73. —3.01*094 217.3 19 187 275 
20200913C 0.18 269.03(6) ^ 65.80(3) 95.5 30.2  576.880(3) 48 43 2 156 + 1.8 (2.85*$12) x 107? 
20200926A 0.22 283.28'010 53.95793 83.7 21.4 758.44(5) 64 59 3 67 £ 43 0.5 (1.58* 128) x 107? 
20180910A 0.43 354.8(9) 89.01(1) 122.6 26.2  684.408(1) 56 55 3 3913 +270 / 4335426 19.1/30.7 (3.83*$23) x 107? 


^ Here we employ the TNS naming convention, see https:/ /www.wis-tns.org/astronotes/astronote/2020- 70. 

> Positions were determined from baseband data (symmetric uncertainties) where available and from per-burst S/N data oth- 
erwise (asymmetric uncertainties; see $3.1), and uncertainties are at 90% confidence. 

* Galactic longitude and latitude for the best position. 

1 Inverse-variance-weighted average DM, and the uncertainty in the calculation of that average. In cases where the DM was 
fixed in the fitburst fit, we set the DM uncertainty to 0.5 pc cm? for the calculation of the average. For candidates with large 
DM variation from burst to burst, this average and its uncertainty may not necessarily be representative and we advise looking 
at the individual burst DMs using generous DM ranges when conducting follow-up searches of these candidates. 

^ Maximum model prediction along this line-of-sight for the NE2001 (Cordes & Lazio 2002) and YMW16 (Yao et al. 2017) 
Galactic electron density distribution models. Neither model accounts for DM contributions from the Galactic halo, which 
contributes no more than 111 pccm ? (Cook et al. submitted). 

f For sources observed twice a day, the second entry corresponds to the less sensitive lower transit. The uncertainties in the total 
exposure for the upper and lower transits of each source are dominated by the corresponding source declination uncertainties 
since the widths of the synthesized beams vary significantly with declination. 

? Fluence completeness limits are given at the 9096 confidence level. For sources observed twice a day, the second entry corre- 
sponds to the less sensitive lower transit. 

^ Adjusted to a 5 Jy ms fluence threshold, assuming a —1.5 power-law energy index. If no burst from a source was detected in 
the time over which the exposure was calculated, no burst rate is reported for that source (84.4). 

! The best-known sky position of this source falls in between the FWHMs at 600 MHz of the synthesized beams. We report the 
exposure at the beam centre of the nearest beam and have scaled the fluence completeness threshold accordingly. 

J Two equally likely “islands” in the localization 9096 confidence region. 


on the main island. In the example of FRB 20190110C (top right), the central region is ruled out by combining 
multiple header localizations, and this is confirmed by the positions derived from baseband data. In the examples of 
FRBs 20191106C (bottom left) and 20191105B (bottom right), combining multiple header localizations reduces the 
allowed phase space in the North-South/declination direction. The baseband positions fall on the mainlobe island for 
the former source, but on a sidelobe island for the latter source. 
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Figure 4. Four representative examples of combined localization patterns, as a function of right ascension (RA) and declination 
(DEC). The gray shaded area shows the 90% confidence region of the combined header (real-time detection) localization with 
the best-known position overlaid in orange. The insets show a zoom-in around the combined baseband localization (orange) 
on top of the individual baseband localizations (purple) shown as well. The combined header localizations allow for multiple 
possible degenerate islands and the baseband localizations can break that degeneracy. 


3.2. Burst and Source Properties 


We characterize each new burst as in the first CHIME/FRB catalog (CHIME/FRB Collaboration et al. 2021). In 
summary, we employ a least-squares fitting algorithm that allows for multiple components and provides measurements 
of DM and scattering (for each burst) and times-of-arrival, widths and spectral properties (for each burst component). 
Width and scattering timescales can be robustly measured from CHIME/FRB intensity data for values larger than 
100 us. Additionally, we determine a bandwidth and duration from the burst dynamic spectrum dedispersed to the 
best-fit DM. Flux and fluence are determined from converting the measured intensities to physical quantities using 
daily observations of bright calibration sources with known flux densities. CHIME/FRB exposure and sensitivity are 
calculated based on metrics collected that included instrument uptime and monitoring of detection S/Ns of known 
pulsars. Finally, we determine a fluence detection threshold for each source to characterize our completeness. This is 
based on a Monte Carlo analysis that takes into account uncertainties in the sky position and the spectral properties 
of the source. Burst rates are calculated from the number of detections and exposure to a source, and linear drift 
rates!? are measured as by Fonseca et al. (2020), except that we do not resample the noise background here — only 
the DM uncertainty. 

Figures 5, 15-17 show dynamic spectra for each burst of each source. As before, most repeater bursts show narrow 
bandwidths of typically 50-200 MHz, with some sources showing bursts with clear cases of downward-drifting subbursts 
(FRBs 20190804E, 20190915D, 20200929C, 20200809E, 20201014B, 20200420A, 20201114A and 20200118D) of a few 
to tens of MHz ms". Interestingly, FRBs 20200127B and 20181226F at modest extragalactic DMs (< 300 pc cm~?) 
show larger bandwidths and narrower widths than is commonly observed for repeaters, reminiscent of bursts from 


12 Using https://github.com/zpleunis/dfdt. 
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other nearby sources FRBs 20181030A (Bhardwaj et al. 2021b) and 20200120E (Bhardwaj et al. 2021a). Some notable 
bursts: FRB 20210327F of source FRB 20201130A shows two resolved components spaced by about 300 MHz and a 
few seconds, FRBs 20200214A, 20200313C and 20200313B of source FRB 20190915D show two subbursts separated 
by ~60, 30 and 30 ms, respectively, and FRB 20201014B of source FRB 20200202A shows at least seven subbursts 
with a spacing that shows a hint of possible periodicity, which is not statistically significant. 

A preliminary long-term periodicity search with an algorithm that uses a Pearson's x? test (e.g., CHIME/FRB 
Collaboration et al. 2020b) for all sources yielded no detections, which is not surprising as the largest number of bursts 
detected for a single source in this sample is only 12. 


3.3. Comparison Analysis 


Whether repeating FRBs represent a distinct population of astrophysical objects from those that have not been seen 
to repeat remains an open question. It is possible that, in time, all FRB sources will eventually emit repeat bursts. 
However, wait times between bursts may be long. It is therefore worthwhile to investigate other possible differences 
between the repeating and apparently-nonrepeating classes of FRBs. 

For this reason, CHIME/FRB Collaboration et al. (2021) reported on a statistical comparison of repeater and 
apparent nonrepeater burst properties. This study revealed significant differences in both burst width and bandwidth, 
as studied in much greater detail by Pleunis et al. (2021b). However, no other statistically significant differences in 
morphological properties (e.g., DM, scattering) were identified. 

Here we repeat the statistical comparison first reported in CHIME/FRB Collaboration et al. (2021), now in- 
cluding the repeater sample reported on here, together with all other published CHIME/FRB-discovered repeaters 
(CHIME/FRB Collaboration et al. 2019a,b; Fonseca et al. 2020; Bhardwaj et al. 2021a; Lanman et al. 2022; Kumar 
et al. 2019). In this analysis, we compare all the same properties as for CHIME/FRB Collaboration et al. (2021), 
between the new, much larger sample of repeating FRB sources, and the Catalog 1 apparent nonrepeaters. We also 
compare the previously published CHIME/FRB repeater sample with the new repeater sample, to compare samples 
discovered with different total sky exposures. The more detailed study of Pleunis et al. (2021b) is beyond our scope, 
and also unlikely very informative as the differences in morphology and spectra clearly persist strongly (see 83.3.2 
below). 

As in the CHIME/FRB Collaboration et al. (2021) analysis, we apply the following constraints when defining our 
comparison samples: 


e Events with bonsai (real-time detection pipeline) S/N « 12 are rejected to mitigate completeness issues. 


e Events with DM < 1.5 max(DMwzaoo1, DMymwie) are rejected to avoid misidentifying rotating radio transients 
or radio pulsars as FRBs due to ~50% errors from models used to estimate the Galactic DM. 


e Events detected in sidelobes are rejected!? since our understanding of the shape of the primary beam at large 
zenith angles is poor, rendering event characterization challenging. 


With these cuts, 155 events are rejected from the nonrepeaters sample, four sources are rejected from the previously 
published CHIME/FRB repeaters sample, and 18 sources are rejected from our new repeaters sample. The remaining 
events represent a total of 305 apparent nonrepeaters, 14 previously published CHIME/FRB repeaters, 26 new repeaters 
(22 from the repeater sample reported here and four’* from other CHIME/FRB-discovered repeaters not included in 
Catalog 1), for a total of 40 CHIME/FRB repeating sources. !* 

Note that below, we also relax some of these constraints in order to probe sensitivity of source selection in our com- 
parison analysis. One significant difference in this analysis compared to that reported by CHIME/FRB Collaboration 
et al. (2021) is that we no longer require the first-detected burst from a repeater to exceed the S/N cutoff, but rather 
that any of the detected bursts exceed the S/N cutoff. This is because our repeaters are now identified through a 
clustering algorithm that can find repeating sources whose first event is not very bright. In our analysis described 
next, we use the first repeating event from each repeating source discovered by CHIME/FRB that satisfies all the 
conditions listed above. 


13 Detections in far sidelobes show multiple narrow knots of emission (~10-MHz wide each, separated by tens of MHz) due to the chromatic 
beam response in the far sidelobes (Lin et al. in prep.). 

1^ FRBs 20190907A (Fonseca et al. 2020), 20200120E (Bhardwaj et al. 2021a), 20171019A (originally discovered by ASKAP; Kumar et al. 
2019) and 20201124A (Lanman et al. 2022). 

15 Tn the Catalog 1 sample, two repeater bursts and four nonrepeaters do not have their flux and fluence measured, as these bursts occured 
right after a system restart with no calibration available and the comparison samples for flux and fluence are thus smaller by six. 
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Figure 5. Dynamic spectra (“waterfall plots"), frequency-averaged time series and time-averaged spectra for all sources, 
ordered from high to low Rec. Bursts from the same source are grouped by the same color. In the upper right-hand corner of 
each panel are the the TNS names and DM values in units of pecm~*. The model fits for which scattering was determined to 
be insignificant are over-plotted in the time series and spectra with thin blue lines, whereas the fits for which scattering was 
determined to be significant are over-plotted in the time series and spectra with thick blue lines. The colored shaded regions in 
the time series and frequency spectra span the respective full width at tenth of maximum (FW'TM) burst widths and emission 
bandwidths. An asterisk in the top right corner of a panel indicates that the event is a potential outlier, based on how it changes 
the Ree of the cluster it belongs to (see 82.2.2). Continued in Figs. 15-17. 
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For distribution comparisons with uncensored data (i.e., those without upper or lower limits), we report both 
Anderson-Darling (AD) and Kolmogorov-Smirnov (KS) tests. A p-value of < 0.01 indicates 99% confidence that two 
samples originate from different distributions. For censored data, we apply statistical survival tests that consider upper 
and lower limits because we lose information if we directly omit censored data. For properties that have upper limits 
(pulse width and scattering time), we use both the Peto & Peto test and the log-rank test, using R NADA package's 
cendiff function!® (Helsel 2005; Lopaka 2020; R Core Team 2020). For events with properties containing lower limits 
(fluence and flux), a right censored survival test is applied using R NADA package’s survdiff function! (Lopaka 
2020; R Core Team 2020; Harrington & Fleming 1982). 


3.3.1. New versus Previously Published CHIME/FRB Repeater Comparison 


Here we compare properties of previously published and new CHIME/FRB repeaters to determine whether the 
sample discovered with longer total exposure time differs from that discovered sooner. For example, a sample discovered 
with longer exposure time might be expected to have a larger mean DM, since a longer exposure may provide more 
distant sources sufficient time to produce two detectable bursts. 

We compare distributions of right ascension, declination, DM, eDM, bonsai S/N, fitburst S/N, boxcar width, 
bandwidth, fluence, flux, width, and scattering. Here, eDM is defined as the DM after subtracting the NE2001 
expected maximum Galactic contribution towards the source. In all cases except flux/?, p-values are > 0.01. However, 
we note that lower limits dominate in the flux values in both the published sample and the new sample (all flux values 
in the published sample are censored). Additionally, the flux values of the new repeaters are more reliable due to better 
localizations while the flux values of the previously published repeaters have more uncertainty. Thus, we cannot be 
confident of the difference in flux. We conclude that there is a lack of significant differences between the two samples' 
properties. This may be due to insufficient difference in exposure time for samples of these sizes to yield a significant 
detection, or too few sources to permit any difference detection. 


3.9.2. Repeater versus Apparent Nonrepeater Comparison 


The p-value results of our statistical tests comparing repeaters and apparent nonrepeaters are shown in Table 2. 
For distributions of right ascension, declination, bonsai S/N, fitburst S/N, fluence, and flux, we find that the 
distributions for repeaters and nonrepeaters are consistent with having been drawn from the same population, as 
we did in our prior analysis. For burst width (both measured and box car) and spectral parameters, we continue 
to see a statistically significant distinction between distributions, consistent with that previously studied by Pleunis 
et al. (2021b), which will be further quantified in a future paper. Pleunis et al. (2021b) found p-value = 10^? when 
comparing width between repeaters and nonrepeaters, and here, we find pap = 1.88 x 107? and pps = 3.73x 1075. For 
bandwidth, Pleunis et al. (2021b) observed p-value = 1075, and we observe pap = 1.24 x 107? and pgs = 5.35 x 1075. 
Pleunis et al. (2021b) obtained p-value — 1074 for boxcar width, and we get pap = 2.49 x 107? and pgs = 6.20 x 107? 
in our analysis. Furthermore, with the new, larger sample of repeaters, we now detect a new potential difference in 
the two distributions, as detailed next. 


DM Comparisons —Figure 6 compares the DM distribution of all one-off events in CHIME/FRB Collaboration et al. 
(2021) and the complete sample of CHIME/FRB-discovered repeating sources of FRBs. Our results suggest a significant 
difference in the DM distributions of repeaters and nonrepeaters. The AD and KS tests give pap = 0.0112 and 
prs = 0.00618, which suggest observed nonrepeaters and repeaters come from different distributions. The analogous 
analysis previously reported by CHIME/FRB Collaboration et al. (2021) using a smaller repeater sample detected no 
such difference. To investigate the putative difference further, we examined how the p-values change as we lower the 
S/N constraint to 11 and 10, in order to increase the sample sizes, though at the possible expense of completeness, but 
noting that the latter applies both to repeaters and nonrepeaters. When the S/N threshold is set to > 11, pap = 0.00425 
and pgs = 0.00324. For S/N > 10, pap = 0.00196 and pgs = 0.00278. Thus, we see that the differences become more 
significant at lower S/N thresholds where the sample sizes increase, as expected for a significant difference. For lower 
S/N cuts, the difference in mean DMs between repeaters and nonrepeaters also increases. At S/N > 12, repeaters have 
a mean DM of 500 + 50pccm~? and nonrepeaters have a mean DM of 662 + 24 pc em ?, with a difference of 162 + 
55 pc cm”? between repeaters and nonrepeaters. At S/N > 11, repeaters have a mean DM of 493 + 47 pccm”? and 


16 https:/ /rdrr.io/cran/NADA /man/cendiff.html 
17 https:/ /rdrr.io/cran/survival/man/survdiff.html 
18 Log-rank test p-value — 0.00357 and Peto & Peto test p-value — 0.0102. 
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Figure 6. Normalized histogram (left) and cumulative distribution (right) of the DM of the 305 nonrepeating events from 


CHIME/FRB Collaboration et al. (2021, blue) and the complete sample of 40 CHIME/FRB-discovered repeating sources of 
FRBs (orange). 
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Figure 7. Normalized histogram (left) and cumulative distribution (right) of the eDM of the 305 nonrepeating events from 


CHIME/FRB Collaboration et al. (2021, blue) and the complete sample of 40 CHIME/FRB-discovered repeating sources of 
FRBs (orange). 


nonrepeaters have a mean DM of 668 + 22 pc cm-^?, resulting in a difference of 175 + 52pccm-^? (means differ with 
> 9996 confidence). At S/N > 10, repeaters have a mean DM of 495 + 45 pc cm ? and nonrepeaters have a mean DM 
of 692 + 23 pc cm ?, giving a difference of 197 + 50 pccm ? (means differ with > 99% confidence). 

This result is also observed in eDM distributions, as shown in Figure 7. With our nominal cuts on the sample, 
we see a possible difference in the underlying distributions, with pap — 0.0105 and pks — 0.00416. The difference 
becomes more significant with lower S/N constraints (for S/N > 11, pap = 0.00464 and pks = 0.00465; for S/N > 
10, pap = 0.00188 and pks = 0.00204), again supporting the existence of a difference in DM distributions between 
repeaters and apparent nonrepeaters. As the S/N threshold is lowered, the difference in mean eDM increases as 
well. For S/N > 12, repeaters have a mean eDM of 436 + 49 pccm ? and nonrepeaters have a mean eDM of 597 + 
24 pccm-^?, with a difference of 161 + 55 pc cm"? between repeaters and nonrepeaters. At S/N > 11, repeaters have 
a mean eDM of 431 + 46 pc cm^? and nonrepeaters have a mean eDM of 602 + 22 pccm-^?, with a difference of 171 


+ 51 pccm”? (means differ with > 99% confidence). At S/N > 10, repeaters have a mean eDM of 431 + 44 pc cm~? 


and nonrepeaters have a mean eDM of 627 + 23 pccm7?, giving a difference of 196 + 50pccm”? (means differ with 
> 99% confidence). 
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Figure 8. Histogram of the scattering times of the 305 nonrepeating events from CHIME/FRB Collaboration et al. (2021, 
light blue) and the minimum scattering times of the complete sample of 40 CHIME/FRB-discovered repeating sources of FRBs 
(light orange) on a log scale. Histograms of scattering times after discarding censored values are overplotted (dark blue for 
nonrepeaters and dark orange for repeaters). 


We thus conclude that with the larger sample of repeaters reported here, we detect a significant difference in DM 
distribution between repeaters and apparent nonrepeaters, with repeater DMs systematically lower than those of 
apparent nonrepeaters. This difference persists when we also include the "silver" sample with the “gold” sample 
(events with Rec < 5) while keeping the same nominal cuts restrictions. The possible reasons for this difference are 
discussed below in 84. 


Scattering Time Comparisons —Next, we look at the distributions of scattering times. As many of our scattering mea- 
surements are upper limits, we make use of statistical survival tests to consider differences in distributions. The Peto 
& Peto and log-rank tests for the repeater versus apparent nonrepeater data indicate similar underlying distributions, 
with p — 0.2 and p — 0.09, respectively. Since large scattering times may be due to confusion with multiple sub-bursts 
merging together, making it hard to identify each sub-burst individually, we also compare the lowest scattering time 
of each repeating source with the scattering time of apparent nonrepeaters (see Fig. 8). We find that the Peto & Peto 
test yields p — 0.1, and the log-rank test yields p — 0.07, respectively, which shows there is no significant difference 
in scattering time distributions in the present samples. Even as we increase the sample size by lowering the S/N 
threshold to 10, we observe no significant difference (i.e., p > 0.01 for the log-rank and the Peto & Peto test). Using 
NADA’s cenfit function!’ (Helsel 2005; Harrington & Fleming 1982; R Core Team 2020; Lopaka 2020), we find the 
Kaplan-Meier estimate mean, which considers upper limits, of the lowest scattering time of repeaters to be 2.32 + 
0.63 ms and that of nonrepeaters to be 3.22 £ 0.50 ms. 


3.4. Burst Rate Analysis 


Here we compute burst repetition rates for all our repeaters, putting them into the broader context of FRB burst 
rates. 

'To compute repetition rates, we first determine the on-sky exposure for each of the sources from the start of the 
experiment, 2018 August 28, through to 2021 May 1. We follow the exposure determination techniques previously used 
by CHIME/FRB Collaboration et al. (2021). In summary, we generate an all-sky exposure map for the aforementioned 
interval by querying metrics which characterize the up-time and sensitivity of the CHIME/FRB system. These metrics 
are combined with the CHIME/FRB beam model, under the assumption that a sky location is detectable if it is within 
the FWHM region of a synthesized beam at 600 MHz. We query the resulting exposure map over a grid of positions 
within the 9096 confidence localization region for each source. The mean and standard deviation of the exposure over 
the corresponding positional grids are provided in Table 1 and Table 3, for the “gold” and "silver" samples, respectively. 
The exposure increases with declination, with the uncertainties in the exposure dominated by the corresponding source 
declination uncertainties. 


19 https:/ /rdrr.io/cran/NADA /man/cenfit.html 


16 THE CHIME/FRB COLLABORATION 


Table 2. Summary of Repeater vs Apparent nonrepeater Distribution Comparisons for Nominal Cutoff Criteria 


Property PAD Dks 
Right Ascension 0.322 0.492 
Declination 0.827 0.914 
DM 0.0112 0.00618 
eDM* 0.0105 0.00416 
bonsai S/N 0.262 0.405 
fitburst S/N 0.691 0.855 
Boxcar width 2.49x 107? 6.20 x 107? 
Bandwidth 1.24 x 1078 5.35x 1078 
Property Peto & Peto | Log-Rank 
Fluence? 0.409 0.217 
Flux? 0.339 0.997 
Width* 1.88 x 1075 3.73 x 107? 
Scattering® 0.222 0.0859 
Lowest Scattering^: 0.140 0.0674 


Notes. 

^ Anderson-Darling probability of originating from same underlying population. 
"Kolmogorov-Smirnov probability of originating from same underlying population. 
Extragalactic DM. 

^Incorporates lower limits — see text. 

“Incorporates upper limits — see text. 

^ Uses only lowest scattering time per repeater. 


In the interval used for the exposure calculation, we detected a total of 128 bursts from the 25 confirmed and 
14 candidate repeaters presented in this paper (i.e., bursts that were detected outside the FWHM at 600 MHz of a 
synthesized beam were excluded from the burst rate computation). While calculating the burst rate, we exclude five of 
these bursts which were detected when the system was not operating nominally. Additionally, for circeumpolar sources 
(6 > +70°), we exclude bursts detected in the lower transit. Therefore, we only compute the burst rate in the upper 
transit for all sources. This choice is motivated by the sensitivity of the primary beam being greater in the upper 
transit, which can provide a stronger constraint on the source activity. The sensitivity to each source also varies based 
on the synthesised beam response. This variation is accounted for in the 9096 confidence completeness thresholds 
reported in Tables 1 and 3. To facilitate comparison across repeaters, we scale all burst rates to a fluence threshold 
of 5 Jy ms, i.e., the average sensitivity of the CHIME/FRB system (Josephy et al. 2021). The scaling involves the 
assumption of a power-law index of —1.5 for the cumulative burst energy distribution of each source. We set the index 
to be —1.5 to be consistent with the repetition rate analyses presented by CHIME/FRB Collaboration et al. (2019b) 
and Fonseca et al. (2020). The scaled repetition rates are plotted in Figure 9. 

We also update the rate estimates for 19 previously published CHIME/FRB repeaters (CHIME/FRB Collaboration 
et al. 2019a,b; Fonseca et al. 2020; Lanman et al. 2022) by computing the exposure up to 2021 May 1 using the 
corresponding source localization regions (see Table 4 in Appendix B). For 13 of these sources, we use unpublished 
baseband localization regions (Michilli et al., in prep.). The localization regions published in the discovery papers are 
used for the rest of the sample. A total of 210 bursts (including unpublished detections) were observed in the upper 
transits of these sources, of which we exclude 15 bursts as they were detected in low-sensitivity periods. We scale 
the repetition rates to a threshold of 5 Jy ms, using the completeness thresholds reported in the discovery papers. To 
allow for comparison of the CHIME/FRB repeaters with the nonrepeaters, we also estimate burst rates for Catalog 
1 one-off sources based on the detection of one burst from each source. Although CHIME/FRB Collaboration et al. 
(2021) report the detection of 474 such sources, we exclude 63 of them because of their detection in low-sensitivity 
periods or in the lower transit. Another 14 sources are excluded as they are found to be repeat bursts from the sources 
presented in this work. The rates for the remaining sources are estimated using the header localization regions and 
completeness thresholds reported by CHIME/FRB Collaboration et al. (2021), and are plotted in Figure 9. 
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Figure 9. Burst rates of previously published and new repeaters, candidate repeaters (the “silver” sample) and Catalog 1 
nonrepeaters as a function of zenith angle (bottom abscissa axis) and declination (top abscissa axis). For each source, the 
uncertainty on the rate is the quadrature sum of the 9096 Poisson uncertainty and the uncertainty in the source exposure (see 
Tables 1 and 3). Histograms of the rate estimates for each class of sources are plotted in the right panel. 


Next we investigate whether the burst rates for the repeating sources are correlated with other properties, namely, 
extragalactic DM, intrinsic width, bandwidth, fluence and scattering. We perform this analysis with the confirmed 
and candidate repeaters presented in this work as well as the 19 previously published repeaters. For these 19 repeating 
sources, we use the burst properties reported in their discovery papers. Additionally, we determine properties of the 
unpublished detections from these sources using the techniques described in Section 3.2. 

While testing for potential correlations, here we define eDM for a source as the mean extragalactic DM estimated 
using the NE2001 and YMW16 models. The source width is assumed to be the average of the intrinsic widths for all 
repeat bursts. If a repeat burst comprises of sub-bursts, we average the widths for the individual sub-bursts before 
using them to calculate the mean width for the source. We obtain the error on the mean width for each source as the 
quadrature sum of the uncertainties in the width measurements of the repeat bursts. We obtain mean bandwidths and 
fluences for each source and the corresponding errors in a similar manner. As for scattering, we adopt the strongest 
constraint among all repeat bursts as the scattering time for each source. This assumes that scattering structures along 
the line of sight to these sources do not change significantly in the interval between detections (but see Ocker et al. 
2022). Due to the finite time resolution of the CHIME/FRB system, some width and scattering time measurements 
are upper limits. In such cases, we assume the 9596 confidence upper limit to be the measured value. 

While checking for a correlation of the repetition rate with these properties, we account for the errors in the 
measurements by performing 10,000 simulation runs. In each run, we sample a possible value of the burst rate and 
also of the property in consideration (say, the width) for each of the repeaters. For each repeater, we sample these 
values from a Gaussian distribution with a standard deviation equal to the lo error on the burst rate and width. In 
each iteration, we check for a correlation using the Spearman's rank-correlation test. We then obtain the median of the 
p-value for the test over all 10,000 iterations. The burst rates and other properties along with the associated median 
p-values are shown in Figure 10. 


4. DISCUSSION 
4.1. Completeness of Repeater Identification 


We believe that our repeater identification is complete for sources with properties similar to previously discovered 
repeaters (e.g., no extreme DM evolution). Next we discuss ways in which our Ree calculation could be made more 
robust and sensitive in the future. 
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Figure 10. Burst rates of previously published and new repeaters (blue) and candidate repeaters (the “silver” sample; green) 
as a function of eDM (top left), mean width (top right), scattering time (bottom left) and mean fluence (bottom right). The title 
of the legend in each panel indicates the associated p-value for the Spearman rank-order correlation (confirmed and candidate 
repeaters; blue and green), with the null-hypothesis being that the plotted quantities are uncorrelated. The p-value for the test 
conducted only using the sample of confirmed repeaters (blue) is indicated in brackets. The errors represent the 90% confidence 
interval for the plotted quantities. 


The clustering analysis (§2.2.1) at the root of the repeater search presented here has used large tolerances in sky 
position and DM, oq = 2.2?/ cos(0), os = 0.5? and opm = 13pccm”?, that reflect uncertainties in the CHIME/FRB 
detection pipeline. It successfully constructs clusters of repeater detections in the main lobes and first sidelobes of the 
beams, but it will miss repeat bursts in the far sidelobes of the telescope, if they are detectable (this depends on the 
luminosity function of the repeaters, see also Lin et al. in prep.). The algorithm groups together events with ADM 
< 13pccm-^?, which makes it sensitive to mild putative DM evolution of repeating sources as well as to a range of DM 
measurement errors in the real-time detection pipeline. We allow for this range because in cases of downward-drifting 
subbursts, the DM can be overestimated by the S/N-optimizing detection pipeline as subbursts get superimposed for 
greater S/N (Hessels et al. 2019). Sources with DM evolution 2 13 pccm ? would have been missed, but this level of 
DM evolution has not been observed in any source of FRBs. 

The Rec calculation is limited by a few observational realities. The probability that we detect an FRB event in 
some volume AaAdADM is a function of CHIME/FRB’s exposure and selection function (82.2.2). This is therefore 
not trivial to calculate. Instead, we empirically estimate it using a normalized histogram of our entire FRB sample, 
including only one burst from each known repeater. This means that the data set that is used to construct the null case 
for the Rec calculation is contaminated by repeating sources. Ultimately, this produces conservative Rec estimates, but 
an improved method would use an uncontaminated null data set or account for the contamination in some way. We 
control for the problem of multiple hypotheses testing by multiplying by N, the total number of FRBs in this sample. 
This à common correction for the look-elsewhere effect, and it has been shown to be conservative (Goeman & Solari 
2014). 
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Another factor that we do not account for is our pipeline-imposed increase in sensitivity for repeaters. Once our 
pipeline identifies a burst as having potentially come from a previously observed FRB, the threshold for saving intensity 
and baseband data in the relevant sky region and DM range is lowered. In principle, this should be accounted for in 
the empirical estimation of Pfa,s,Dm) but is non-trivial and will be the subject of future work. 

The time between detections is also not accounted for in our Ree calculation. Intuitively, FRBs with similar DMs 
and spatially coincident localizations that are detected in rapid succession (e.g., within a single CHIME/FRB transit) 
may seem more likely to be repeaters than those detected years apart. However, how to include this timing is not 
clear, particularly because we do not know the general temporal behavior of repeaters. More detailed quantitative 
characterizations of the waiting time distributions and other temporal properties of the repeater population could be 
helpful in this regard. 

Future improvements to the Rec framework could also allow for priors on additional parameters, such as the burst 
duration and bandwidths, which have been shown to be different for repeater bursts and one-off events (Pleunis et al. 
2021b), or the scattering properties, when we understand better what variability can exist (see, e.g., Ocker et al. 2022). 
Regardless, CHIME/FRB aims to provide sky positions with subarcminute precision for many future FRBs through 
its “Outriggers” upgrade (Leung et al. 2021; Cassanelli et al. 2022), which will make identifying true repeaters among 
the sample easier in the future. 


4.2. Comparison of New and Previously Known Repeater Samples 


As stated in 83.3.1, we conclude that there is no significant difference between new and previously published 
CHIME/FRB repeater samples. All p-values are > 0.01 except flux, which has a log-rank test p-value « 0.01 but has 
large uncertainty due to the presence of many censored values. We might expect to detect more repeaters in this new 
sample at lower declination where CHIME/FRB’s sensitivity is lower, since more observation time has passed, but no 
difference in the declination distributions is observed. Similarly, due to longer total exposure, we might expect to see 
a difference in the DM distributions between previously published and new samples of repeaters since there is more 
opportunity to detect more distant repeaters, but no difference is observed. 

This lack of any apparent difference could indicate no difference in the two samples’ properties, insufficient difference 
in exposure time for samples of these sizes to yield a significant detection, or too few sources to permit detection of any 
differences. Alternatively, the lack of an apparent difference between the new and previously published repeaters could 
be a result of FRBs having finite lifetimes, since the active lifetime of an FRB could end before a burst is detected. 
Or, new repeating sources (cf. FRBs 20201124A and 20220912) could “turn on” at all declinations and DMs at a rate 
comparable to the increase in discovery of sources at high DM and low declination with increased total exposure. The 
two competing rates could mask the expected increase in detections of more distant sources over time, especially when 
the total number of sources is still relatively low. 


4.3. Comparison of Apparent Nonrepeating and Repeating Sources 


As discussed in §3.3.2 and shown in Table 2, the AD and KS tests indicate that the DM and extragalactic DM 
(eDM) distributions of the full repeaters sample and Catalog 1 apparent nonrepeaters (meeting the criteria specified 
in §3.3.2) are seemingly inconsistent with arising from the same population. This inconsistency increases as we lower 
the S/N threshold. A likely explanation for the latter observation is that the sample size increases as we lower the 
S/N threshold, making any real differences more significant. If so, the difference in DM distribution between repeaters 
and apparent nonrepeaters is very significant. Such a difference could in principle indicate new, important support for 
the two source classes being of physically different origins. A difference in DM distribution could indicate a difference 
in one or more of (i) distribution in redshift space, (ii) host galaxy type, or (iii) local environment. However, none of 
these is necessary to explain the effect. Gardenier et al. (2021) showed that the detection of nonrepeaters and repeaters 
is inherently subject to different observational selection effects. If repeaters and nonrepeaters originate from the same 
population, we would observe a difference in DM distributions because a second bright burst required for a source to 
be identified as a repeater is more challenging to detect at higher DM, if the luminosity function universally declines 
at higher luminosity. Detailed population modelling will be required to determine whether the degree of the difference 
reported here can be explained by this simple bias assertion alone under reasonable assumptions of luminosity function, 
or whether one can conclude there is a genuine difference in source classes. 

We do not observe a difference in scattering time distributions when we compare the first detected burst per source 
after the cuts stated in §3.3.2, nor do we observe a difference between the two distributions when we compare the 


20 THE CHIME/FRB COLLABORATION 


lowest scattering time of each repeating source with nonrepeaters. We note that it is possible for scattering-time 
measurements to be biased by broad pulse widths and complex morphology. Since repeater bursts have been shown to 
have longer duration and more complex morphology (Pleunis et al. 2021b), we might expect a bias in which repeater 
bursts have longer scattering times compared to nonrepeaters. Yet, we observe repeaters to have comparable or even 
possibly lower scattering times both when we compare the first burst that passes the criteria in 83.3.2 (repeaters: 2.57 
+ 0.77 ms; nonrepeaters: 3.22 + 0.50 ms) and when we compare the minimum scattering time of each source (repeaters: 
2.32 + 0.63; nonrepeaters: 3.22 + 0.50ms). On the other hand, if there were a strong DM/scattering correlation, 
then repeaters having a lower mean DM and eDM than apparent nonrepeaters, as we observe, would suggest they 
should also have lower scattering times. However, strong DM/scattering correlations have not been seen in the FRB 
population (e.g., Chawla et al. 2022). Detailed population studies, preferably on burst properties measured using 
voltage data, for which the time resolution is higher, hence scattering may be more easily measured, will be needed to 
properly search our larger sample for such a correlation. 


4.4. Burst Rates 


Here we discuss the burst rate analysis shown in 83.4, where we found burst rates for all new (both *gold" and 
“silver” samples) and previously published CHIME/FRB repeaters, as well as for Catalog 1 nonrepeating sources (see 
Fig. 9). 

First, we note that the estimated rates might not be equal to the true repetition rates of these sources because of 
two assumptions made in this analysis. First, assuming a power-law index of —1.5 for scaling the burst rates to a 
common fluence threshold will introduce rate errors of varying magnitudes for different sources. This is because the 
measured index for known repeaters varies widely, ranging from —0.7 to —3.6 (Law et al. 2017; Lanman et al. 2022). 
Second, some of the bursts in our sample could have been detected in the sidelobes of the synthesized beams. As 
the exposure map does not include exposure in the sidelobes, the true repetition rate for some sources could be lower 
than is reported here. However, as both these assumptions affect the estimated rates across all sources in a stochastic 
manner, they are unlikely to affect the conclusions of a rate comparison between the repeating and non-repeating 
populations and among the repeating sources. 

Among the repeating sources, we find that the candidate repeaters have, on average, lower repetition rates as only 
two bursts have been detected from most of these sources. We also find that the burst rate is, on average, lower for 
high-declination sources. This is because the increased exposure at high declinations allows us to detect repeating 
FRBs with low repetition rates and place stricter constraints on the rates for the one-off sources. While some active 
repeating sources, such as FRB 20201124A (Lanman et al. 2022), have anomalously high repetition rates, we find no 
clear bimodality in the rates between repeating and apparently non-repeating sources (right panel of Fig. 9). Therefore, 
we cannot as yet rule out the possibility that all FRBs repeat. This conclusion is strengthened by the observation that 
some of the candidate repeaters have lower rates than the rate upper limits of one-off sources. The transit path for 
these candidate repeaters has higher sensitivity than the all-sky average, which results in the scaled detection rate (for 
a 5-Jy ms threshold) being lower than limits for one-off sources. The fact that we are preferentially seeing repeaters at 
sky locations with the highest sensitivity and exposure suggests that all FRBs could be seen to repeat, given enough 
exposure and sensitivity. It is beyond the scope of this paper to investigate how parameterizing clustered repetition 
would affect the results of this analysis, but would be an interesting future investigation. 

In 83.4, we also searched for correlations between burst rate and various burst properties (see Fig. 10). The burst 
rate might be expected to be anti-correlated with eDM if the repeater luminosity function is, on average, steep, and if 
eDM is on average higher for more distance sources, both likely possibilities. Higher rates for lower eDM under these 
assumptions would be expected because closer sources and their repeat bursts will be easier to detect. Although we 
see a slight hint of this for the confirmed repeaters (p-value — 0.05 for confirmed sources, but only 0.2 for the full 
sample; see Fig. 10 top left), overall, the correlation is weak, suggesting more statistics may be needed to detect such 
an effect with high significance. 

On the other hand, we see a significant correlation between average pulse width (defined as the average width of all 
detected bursts for each source) and burst rate (p-value — 0.00005 for the full sample, and 0.02 for just the confirmed 
sources; see Fig. 10 top right). The CHIME/FRB instrumental selection function strongly disfavors detection of long 
bursts (CHIME/FRB Collaboration et al. 2021; Merryfield et al. 2022), so preferentially detecting more wide ones 
cannot be due to that bias. On the other hand, CHIME/FRB’s detection pipeline operates on data having time 
resolution of 0.983 ms, so much narrower bursts, especially if faint, are selected against, a bias not made clear by the 
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constant fluence injection model used when demonstrating our bias (see Merryfield et al. 2022). Hence, bursters having 
very narrow, frequent bursts (e.g., Nimmo et al. 2022b) are not well represented here. Nevertheless, the correlation we 
detect is present in our data and may represent a genuine astrophysical trend. Such a correlation might be expected 
as longer-duration bursts may extend over wider solid angles, hence be more detectable (e.g., Connor 2019). 

We also detect a possible positive correlation between burst rate and scattering time (here defined as the smallest 
scattering time among the repeats from each source) with p-value — 0.004 for the full sample, but only 0.1 for the 
confirmed sources (see Fig. 10 bottom left). As rate is likely an intrinsic source property, whereas scattering time 
is likely related to the propagation path, a true correlation would be somewhat surprising if scattering originates 
anywhere but the environment local to the source. On the other hand, such a correlation could arise if, for example, 
a more active source is preferentially embedded in a denser, more turbulent plasma. We note that highly scattered 
bursts are strongly selected against in our CHIME/FRB detection pipeline (CHIME/FRB Collaboration et al. 2021; 
Merryfield et al. 2022), so detecting more at large scattering time is unexpected given this instrumental bias. On the 
other hand, many of our scattering times are upper limits (here treated as measurements), either because the time 
is smaller than our time resolution, or because it is smaller than the burst width, so the associated p-values must be 
treated as approximate at best. Nevertheless, a possible correlation is present in our data but seems likely related 
to the rate/width correlation discussed above, since, again, detection of scattering is often strongly limited by pulse 
width. We therefore do not ascribe great significance yet to this possible correlation and resist any astrophysical 
interpretation for the time being. 

Finally, we detect a hint of positive correlation between burst rate and mean burst fluence (p-value 0.006 for confirmed 
sources and 0.07 for the full sample; see Fig. 10 bottom right). This could be a selection bias; for example, closer 
sources will tend to have higher fluences, and hence be easier to detect repeatedly. On the other hand, the absence of 
as strong a correlation between burst rate and eDM does not support that hypothesis. Alternatively, there could be a 
correlation between burst rate and luminosity, wherein the most active sources are also the most radiatively energetic. 
On the other hand, there could be no such correlation, yet we tend to detect the most luminous, active sources first. 
If so, the top left corner of the plot will eventually fill out. 

In summary, a preliminary study of correlations between repeater burst rate and various burst properties reveals 
some possible relationships that are intriguing. We note that in some cases, the correlations are significant with just the 
confirmed sample and sometimes with the full sample (i.e., including the “silver” sample presented here), suggesting 
that one or more of the putative correlations may be spurious. CHIME/FRB will continue to detect repeating sources 
(see 84.5 below), which will help improve the statistics. Detailed simulations accounting for the variety of instrumental 
and observational biases at play will be required before astrophysical conclusions can be drawn. However, these are 
underway. 


4.5. Repeater Fraction 


If all FRBs eventually repeat, the number of repeaters over the total number of observed FRB sources, frep, will 
tend to an equilibrium close to unity as FRB surveys continue. The value of the equilibrium depends on the birth rate 
and lifetime of repeaters. If only some FRBs repeat, frep will plateau or turn over at a value significantly different 
from 1 (Ai et al. 2021; Gardenier et al. 2021). The evolution of fep furthermore depends on the luminosity function 
and repetition rate distribution of sources. To study frep, we split our experiment into five declination bins with 
approximately equal exposure and sensitivity within each bin. We separate upper and lower transits around the 
NCP (bursts with circles and triangles as markers, respectively, in Fig. 3) which effectively results in six independent 
declination bins. T'his split is necessary because we want to investigate the dependence of the repeater detection rate 
on exposure and sensitivity; with a uniform distribution of repeaters on the sky, it is possible that we run out of new 
repeaters to detect more quickly at sky positions with more daily exposure. 

To ensure completeness, we only include events in this analysis with detection S/N > S/Ntnreshola, with S/Nthreshola € 
{10,11,12}. For this analysis, we “discover” a repeating source of FRBs the second time we detect from it a burst 
with S/N > S/N esnoia (highlighted with an extra circle or triangle in Fig. 3). We count the number of FRBs, Nprp, 
and repeating source discoveries, Nrep, up to each UTC day and assign Poisson counting errors (ow 2 VN, for large 
N, using scipy.stats.poisson.interval to be precise, also for small N). We extrapolate the FRB detection rate 
after the cutoff date for the first catalog (2019 June 30) by approximating the detection rate as the total number of 
FRBs detected in each analysis bin divided by the exposure up until the cutoff date. 
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from propagating errors. For the uncertainty calculations, we keep Nprp and on,,, fixed to their measurements on 
2019 June 30 instead of calculating uncertainties based on extrapolated values. We report the 9096 confidence interval. 

The repeater fractions are shown as a function of average exposure (in a given declination bin) in Figure 11 for a 
S/N > 12 threshold and without imposing a cutoff of the maximum extragalactic DM. It can be seen that frep tends 
to an equilibrium of 1-4% in all bins, with consistent values between bins at the 90% confidence level. Summing 
FRBs and sources over all bins, we find an overall frep = 2.6153%. The repeater fraction initially reached ~20% 
in the 49.6? < 6 « 69.8? bin because FRB 20180916B (CHIME/FRB Collaboration et al. 2019a) was discovered 
as a repeating source early on. The computing node that processes eight beams in the 27? « 6 « 33° range has 
been sporadically on and off during periods of the experiment. We have thus checked that excluding this declination 
range does not affect our results significantly, and have eventually included this declination range for the measurement 
presented here. We do not currently have the statistical power to differentiate between frep tending to an equilibrium 
or turning over. We have also inspected the FRB and repeater detection rate as a function of UTC time (not exposure, 
as in Fig. 11) and find no evidence for a change of rate there either. The rates as a function of exposure and as a 
function of UTC time probe slightly different properties of the sources, as the total exposure in all cases is collected 
over a ~three-year period. 

We investigate the behavior of frep as we lower the S/N threshold and as we probe only a local volume of the 
Universe by imposing a maximum extragalactic DM for sources to include. The evolution is shown in Figure 12. At 
maximum extragalactic DM < 200pccm-?, we have not yet discovered any repeating sources in some bins, likely 
due to small number statistics. In general, we observe that frep is higher at lower maximum extragalactic DM and 
gradually falls off as we probe higher extragalactic DMs (in line with findings in 83.3.2 that repeaters have lower mean 
DM and extragalactic DM). However, in all cases, the repeater fractions in different declination bins are consistent 
at the 90% confidence level. We observe a higher frep in the two bins around CHIME's zenith (ô =49.6°) where our 
sensitivity is highest. This could just be a selection effect from being more sensitive to fainter repetitions and thus 
detecting repeaters at a higher rate. The trends are similar regardless of S/N threshold and the use of the NE2001 or 
the YMW16 model for the estimation of the extragalactic DM. 

Shin et al. (2022) used a detailed population synthesis to infer a volumetric rate of of 7.3 x 10“ bursts Gpc^? year” 
above a pivot energy of 10?? erg and below a scattering timescale of 10 ms at 600 MHz based on FRBs in CHIME/FRB's 
first catalog. Based on the observed frep in CHIME/FRB, the rate of active repeating sources can only be a few percent 
of that rate. Without an analogous population synthesis accounting for the observed frep, we cannot make any firm 
conclusions about the statistics of the population of repeating FRBs versus a putative distinct population of truly one- 
off events. It is also not possible to compare frep as measured for the CHIME/FRB survey with frep in other surveys 
because other repeating sources of FRBs were discovered in targeted follow-up observations rather than untargeted 
observations (Spitler et al. 2016; Kumar et al. 2019; Luo et al. 2020; Kumar et al. 2021; Niu et al. 2022). 


1 


4.6. Do All FRBs Repeat? 


Even though multiple lines of evidence point to observed differences between the properties of repeater bursts and 
apparently one-off events, it is not yet possible to rule out that all FRBs repeat. Detailed population synthesis studies 
are necessary to properly interpret observations and rule out selection effects as the origins of the observed differences. 

'There is an observed dichotomy between the burst durations and bandwidths of repeater bursts and apparently 
one-off events (Pleunis et al. 2021b) that persists in this new sample (83.3). To this, we add an observed difference 
in DM distributions (83.3.2) but no detectable bimodality in repeat rates (84.4). We find a significant correlation 
between burst width and repetition rate only if we include candidate repaters in the analysis, which might point to an 
evolution of burst properties with rate (84.4). Only 2.6% of sources detected by CHIME/FRB has so far been seen to 
repeat. 
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69.8? <6 < 90.0? (L) 


69.8? « 5 « 90.0? (U) 


49.6? < 6 < 69.8? 


29.4? « ó « 49.6? 


9.2? « ó « 29.4? 


—11.0? < ô < 9.2? 


0 2 4 6 8 10 


Average exposure (days) 


Figure 11. Fraction of repeating FRBs over all FRB sources, frep, as a function of average total exposure in different declination 
bins (ranges are indicated in the top left or right corners, with U/L separating the “upper” and “lower” transit of the top bin, 
respectively) with approximately equal exposure and sensitivity within each bin. In each pair of panels the first shows the 
number of FRB (green) and repeater (orange) sources and the second value of frep (purple) with the values at 2021 May 1 
written out (90% confidence). The FRB rate is extrapolated after the cutoff date (gray dotted line) for the first catalog. This 
figure includes only detections with S/N > 12. 
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Figure 12. Evolution of frep (in 96) in six declination bins as a function of maximum extragalactic DM (NE2001) for three 
different S/N thresholds. While there are hints of trends towards higher frep at lower maximum extragalactic DM and closer to 
zenith (ô =49.6°), note that uncertainties are not shown here and that all values are consistent at the 90% confidence level. 


Reconciling the dichotomy in burst morphology with one population of repeating FRBs with a continuum of repetition 
rates needs a correlation between repetition rate and burst duration (potentially observed; 84.4) and anti-correlation 
between repetition rate and bandwidth (not yet observed). This could either be achieved intrinsically through the 
emission mechanism or extrinsically through a propagation effect (e.g., tied to an evolutionary stage if active repeaters 
are young sources). Beaming geometry could provide a correlation between repetition and duration as wider opening 
angles could make bursts from sources both easier to detect and longer in duration (Connor et al. 2020). However, it 
is unclear how a wider opening angle could also lead to narrower emission bandwidths. 

Studies of the host galaxies (Bhandari et al. 2022) and galaxy offsets (Mannings et al. 2021) of FRBs have not 
found clear differences between those of repeating sources and one-off events, except that the hosts of repeaters 
FRBs 20121102A and 20190520B are outliers in the currently known population. 

An additional future probe of differences could be polarization properties, as more polarization properties of one- 
off events become available. On the other hand, there already exists an observed wide diversity in the polarization 
properties of the repeaters (e.g., Mckinven et al. 2022b). 

Matching the observables to an astrophysical population of FRBs requires population synthesis. This has been used 
previously to constrain the repetition statistics and energy function of repeaters (Caleb et al. 2019; James et al. 2020a,b; 
Lu et al. 2020) and, for CHIME/FRB data, to interpret the scattering timescales of FRBs in the first catalog (Chawla 
et al. 2022). A publicly available general framework for population synthesis is provided by frbpoppy”? (Gardenier 
et al. 2021). Here, we outline some considerations for accurately modeling FRBs with CHIME/FRB observations for 
future population studies. 

When modeling repeaters and one-off events as distinct populations, many new variables are introduced. Apart from 
separate volumetric rates, it needs to be carefully considered whether they share other population model parameters, 
such as those describing their energy function, evolution through cosmic time, and spectral properties. Shin et al. 
(2022) fit an energy function and volumetric rate to CHIME/FRB’s first catalog. As ~90% of FRBs in the first catalog 
are one-off events, it is likely that those measurements could be used well to model one-off events. Those and other 
model fits are, however, sensitive to the inclusion of repeating sources for which only one event has been detected 
as apparent nonrepeaters. There is also the question of how to model the burst rates of repeaters, as varying levels 
of burst activity have been observed across the population of identified repeating FRBs, including highly clustered 
and rarely repeating sources (§4.4; see also, e.g., Oppermann et al. 2018; Good et al. 2022; Nimmo et al. 2022b). 
Population synthesis will need to reproduce the observed distribution of repetition rates, while also reproducing the 
observed differences in other burst properties. 


20 https: //github.com/davidgardenier/frbpoppy 
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In addition to modeling the underlying intrinsic population, consideration also needs to be taken towards modeling 
the survey in which these bursts were detected and comparing the modeled population bursts with the observed data. 
For this, modeling selection effects is crucial. To accurately forward-model the CHIME/FRB survey, it is important 
to take into account its exposure?! and to use a beam model.? The CHIME/FRB synthesized beams have chromatic 
sensitivities set by the primary beam response of the CHIME telescope and the beamforming algorithm, which can lead 
to, e.g., apparent, observed, bandwidths that are smaller than intrinsic bandwidths if bursts are detected away from 
the centers of beams (see, e.g., 82.1.3 and Fig. 2 in Pleunis et al. 2021b). The selection function of the CHIME/FRB 
detection pipeline has been characterized with injections and has been made publicly available (Merryfield et al. 
2022).2° The strongest bias of the detection pipeline is against bursts with scattering timescales above 10ms at 
600 MHz (CHIME/FRB Collaboration et al. 2021). 


5. CONCLUSIONS 


We have presented 25 new repeating sources of FRBs, plus 14 additional candidate repeaters, from a complete search 
of the CHIME/FRB database from the start of operations on 2019 September 3 up to 2021 May 1. This brings the 
total number of repeating FRBs discovered by the CHIME/FRB project to 46, or counting the candidates, 60. This 
large bounty is a result of CHIME's large field-of-view, high sensitivity, and daily surveying of a large area of sky. 
Additional repeater discoveries continue to be made by CHIME/FRB at a steady pace. 

The DMs of the confirmed sources reported on in this paper fall in the range 200-1700 pc cm^?. Combining these 
with the DMs of previously reported CHIME/FRB repeaters, we have now demonstrated that on average, the DMs 
of repeaters are lower than those of the apparently nonrepeating sources found by CHIME/FRB. This seems likely 
to be at least partly an observational bias, since, for reasonable luminosity functions, nearby repeaters will be found 
preferentially. However, it may also signal a distinct astrophysical origin for repeaters versus nonrepeating sources; 
detailed population synthesis studies that account for the many observational and instrumental biases will be required 
to disentangle such effects. The possibility of distinct physical origins is supported but unproven by statistically 
different dynamic spectra between the two groups — wider bursts and narrower spectra for repeaters — as previously 
shown (CHIME/FRB Collaboration et al. 2021; Pleunis et al. 2021b), a result which is strongly supported by our 
analysis of our now much larger sample. 

Interestingly, we have detected no significant bimodality in the repeat rates of repeaters and apparent nonrepeaters, 
although the most active CHIME/FRB repeaters do have anomalously high rates compared to the bulk of the studied 
population. While we measure that 2.613% of sources detected by CHIME/FRB repeat, the lack of bimodality in 
burst rates between the populations suggests that some fraction of our apparent nonrepeaters may yet repeat and that 
we cannot rule out all FRBs eventually repeating. 

We encourage radio and multi-wavelength follow-up of the newly presented sources, and especially of our candidates. 
As for previously published repeaters, we update a webpage?^ with new detections of these sources to communicate 
which sources are active while also continuing to provide the community with realtime events via VOEvents.?? 


21 The exposure up to 2021 May 1 is available here: https://www.canfar.net/storage/vault/list/ AstroDataCitationDOI/CISTI. CANFAR/23. 
0004/data/exposure 


22 An early version of our beam model, used for injections, is publicly available here:  https://github.com/chime-frb-open-data/ 
chime-frb-beam-model 


23 https:/ /www.canfar.net/citation/landing?doi—22.0005 
24 https:/ /www.chime-frb.ca/repeaters 


25 https:/ /www.chime-frb.ca/voevents 
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APPENDIX 


A. CANDIDATE REPEATING SOURCES OF FAST RADIO BURSTS 


Source properties of 14 candidate repeating source of FRBs are summarized in Table 3. The dynamic spectra of 
their bursts are shown in Figure 13. These candidates have lower significance, which by definition means the burst-to- 
burst DM and sky position differences are larger and the repetition rates are low. The DM differences are as large as 
~14pcecm ^? (candidate FRBs 20190107B and 20200320A) and variation of this magnitude has not yet been observed 
in repeating sources of FRBs before. Some of the bursts in this candidate sample show larger than average bandwidths 
for repeater bursts, which we caution to overinterpret (i.e., we do not downweight their likelihood as candidates for this 
reason and we do not claim a correlation between bandwidth and repetition rate using this candidate sample). Some 
candidate sources also show apparently different scattering timescales from burst to burst (candidate FRBs 20190328C, 
20201105A, 20200828A, 20210203E and 20190127B) that we cannot at this time confirm are real. Scattering variations 
have been detected in FRB 20190520B (Ocker et al. 2022). 


Table 3. Properties of 14 New Candidate Repeating Sources of FRBs, Ordered by Increasing Rec (Our “Silver” Sample) 


FRB Source” Ree a? ó^ Ie b° DM? DMxzg»o: DMYuwie Noursts Exposure! Completeness? Burst rate” 
(J2000) (J2000) (deg) (deg) (pccm^?) (pc cm~?) (pe cm~’) (hr, upper / lower) (Jy ms) (hr?!) 
20190303D 0.62 185.34(7)  70.69(2) 126.5 46.2 714.552(8 37 29 2 55-£44/142-£3 2.8 / 16.7 (1.49 22) x 107? 
20181201D 0.64  267(1) 89.12(1) 122.0 27.4 443.849(2 54 5 2 3591 + 288 / 4359 + 14 2.7 / 4.2 (2.24* 181) x 1074 
20190328C | 0.71 75.657019  82.11*072 130.6 23.3  472.819(8 64 67 2 3114161 / 294 + 161 1.8 / 5.7 (7.16*27:10) x 1074 
20191105B! 0.73 — 8411(6) 72.52(2) 140.6 20.4  311.226(2 74 84 2 1633-15 / 161 £0 3.1/ 511 (5.9*15;7) x 107? 
20190107B/ 0.78  49.31*]!5 83.407013 127.1 21.8  168.253(3 68 73 2 362 + 189 / 329+ 196 0.2 / 0.4 (1.91* 220) x 1075 
22.37 02, — 83.377029 124.1 20.6 72 8 

20200320A 0.92 42.45*002  15.84*007 160.2 —38.3  593.524(2 46 39 2 5911 8.7 (3.91* 1458) x 107? 
20190210C 1  313.90'012 89.107021 122.2 26.7 643.365(2 54 5 2 3354 + 1180 / 3086 + 959 1.7 / 2.6 1.83x107* 
20190812A 1.4 268.0502? 50.71*002 78.1 29.8 252.889(2 48 4 2 81 +26 0.9 (1.85* 103) x 107? 
20200828A 1.6 114.08*024 80.027009 134.1 28.6 — 555.41(2) 54 50 2 241+ 133 / 223 + 129 0.4 / 1.6 (1.89* 128) x 1074 
20190905A — 2.2 73.27499  88.24*'007 124.6 26.2  231.5577(7) 57 55 2 5399 + 980 / 3928 + 1313 0.4 / 0.6 (9.9922?) x 1079 
20190127B 2.7  169.26*006  83.51*005 126.0 33.0 678.350(3 48 42 2 345 + 191 / 358+ 197 0.1 / 0.3 (7127) x 1078 
20210323C 3  12207*011 72.35*039 142.6 31.6 288.6(3) 50 45 3 135 78 / 108 + 60 3.6 / 31.8 (1.373239) x 107? 
20180909A 3.9 120.044932 57.007009 160.6 317 418.52(4) 52 49 2 75 +45 0.1 (3.35* 142) x 1075 
20200508H 4 54.600229  88.04*022 1244 25.7 479.430(6 58 57 2 1296 + 537 / 1395 + 375 0.8/11 (1.04*227) x 107“ 


^ Here we employ the TNS naming convention, see https:/ /www.wis-tns.org/astronotes/astronote/2020- 70. 

> Positions were determined from baseband data where available and from per-burst S/N data otherwise (see 83.1), uncertainties 
are at 9096 confidence. 

^ Galactic longitude and latitude for the best position. 

1 Inverse-variance-weighted average DM, and the uncertainty in the calculation of that average. In cases where the DM was 
fixed in the fitburst fit, we set the DM uncertainty to 0.5 pc cm"? for the calculation of the average. For candidates with large 
DM variation from burst to burst, this average and its uncertainty may not necessarily be representative and we advise looking 
at the individual burst DMs using generous DM ranges when conducting follow-up searches of these candidates. 

^ Maximum model prediction along this line-of-sight for the NE2001 (Cordes & Lazio 2002) and YMW16 (Yao et al. 2017) 
Galactic electron density distribution models. Neither model accounts for DM contributions from the Galactic halo, which 
contributes no more than 111 pccm ? (Cook et al. submitted). 

f For sources observed twice a day, the second entry corresponds to the less sensitive lower transit. The uncertainties in the total 
exposure for the upper and lower transits of each source are dominated by the corresponding source declination uncertainties 
since the widths of the synthesized beams vary significantly with declination. 

? Fluence completeness limits are given at the 9096 confidence level. For sources observed twice a day, the second entry 
corresponds to the less sensitive lower transit. 

^ Adjusted to a 5 Jy ms fluence threshold, assuming a —1.5 power-law energy index. If no burst from a source was detected in 
the time over which the exposure was calculated, no burst rate is reported for that source (84.4). 

* The best-known sky position of this source falls in between the FWHMs at 600 MHz of the synthesized beams. We report the 
exposure at the beam centre of the nearest beam and have scaled the fluence completeness threshold accordingly. 

J Two equally likely “islands” in the localization 90% confidence region. 
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B. BURST RATES FOR PREVIOUSLY PUBLISHED REPEATERS 


'Table 4. Updated Burst Rates for Previously Published Repeaters 


FRB Source Reference Burst Rate^ 
(hr?) 
20180814A — CHIME/FRB Collaboration et al. (2019a) —(6.57*$92) x 107! 
20180916B ^ CHIME/FRB Collaboration et al. (2019b) — (4.48* 590) x 107? 
20181030A — CHIME/FRB Collaboration et al. (2019b) —(1.48*727) x 107? 
20181128A — CHIME/FRB Collaboration et al. (2019b) —(4.08*2$$8) x 107? 
20181119A — CHIME/FRB Collaboration et al. (2019b) —(6.77*$29) x 107? 
20190116A = CHIME/FRB Collaboration et al. (2019b) —(3.18*$83) x 107? 
20181017A — CHIME/FRB Collaboration et al. (2019b) - 
20190200A — CHIME/FRB Collaboration et al. (2019b) (6.95*15,:29) x 107? 
20190222A — CHIME/FRB Collaboration et al. (2019b) — (1.44*299) x 107? 
20190208A Fonseca et al. (2020) (43475 5 x 107? 
20190604A Fonseca et al. (2020) (7.961110) x 1073 
20190212A Fonseca et al. (2020) (3.364398) x 107? 
20180908B Fonseca et al. (2020) (diit x 107? 
20190117A Fonseca et al. (2020) (2.1572:39) < 107! 
20190303A Fonseca et al. (2020) (124679 1. x 107! 
20190417A Fonseca et al. (2020) (10975 3) x 107? 
20190213A Fonseca et al. (2020) (3847558) x 107? 
20190907A Fonseca et al. (2020) (2.98*293) x 107? 
20200120E Bhardwaj et al. (2021a) - 
20201124A Lanman et al. (2022) 3.281166 
20171019A Kumar et al. (2019) - 


^ Adjusted to a 5 Jy ms fluence threshold, assuming a —1.5 power-law energy index. If no burst from a source was detected in 
the time over which the exposure was calculated, no burst rate is reported for that source (84.4). 
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Figure 13. Same as Figure 5, but for candidate repeating sources of FRBs (i.e., “silver” sample). 
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C. CATALOG EXCERPT 


In Table 5, we provide an excerpt from our catalog of all repeater candidates presented in this work. We show all 
fields for four FRB entries. Field descriptions can be found in Table 6 and the full data set accompanies the online 
version of this article and at the CHIME/FRB Public Webpage.”° 


Table 5. Excerpt from Repeater Catalog 


tns name previous name repeater name rai ra | err low ra i err up 
(degrees) (degrees) (degrees) 
FRB20191106C v 199.577 0.013 0.013 
FRB20200621D FRB20180910A 354.83 0.87 0.87 
FRB20190804B FRB20190303D 185.339 0.067 0.067 
FRB20210224A FRB20201221B 124.199 0.026 0.026 
ra_2 ra_2_err_low ra_2_err_up ra_notes dec_i dec_i_err_low 
(degrees) (degrees) (degrees) (degrees) (degrees) 
43.0025 0.0088 
89.014 0.011 
70.691 0.024 
48.781 0.017 
dec. 1 err up dec.2 dec.2.err low dec.2. err up dec.notes gl 
(degrees) (degrees) (degrees) (degrees) (degrees) 
0.0088 105.7 
0.011 122.59 
0.024 126.52 
0.017 170.55 
gb exp_up exp.up.err exp.up.notes exp.low exp.low err 
(degrees) (hours) (hours) (hours) (hours) 
13.22 65.6 2.9 "m ia 
26.19 3910 270 4335 26 
46.23 55 44 141.7 2.7 
33.81 90.3 1.6 m T EE 
exp.low notes bonsai_snr bonsai_dm low_ft_90 up_ft_90 low_ft_95 
(pe em?) (Jy ms) (Jy ms) (Jy ms) 
9.36 331.6 0.0 1.9 0.0 
29.97 695.5 30.7 19.1 38.0 
9.33 716.6 16.7 2.8 21.2 
E 9.22 516.8 0.0 3.4 0.0 
up. ft 95 snr fitb dm fitb dm fitb err dm. exc.1.ne2001 dm exc.2.ne2001 
(Jy ms) (pc cm?) (pc cm 3) (pc cm?) (pc cm?) 
2.4 19.83 333.057 0.027 308.0 
24.3 23.2 700 700 640.4 
3.4 10.43 716.0068 0.0091 679.4 
4.6 20.95 510 510 458.9 ahs 
dm_exc_1_ymw16 dm exc.2 ymw16 bc.width scat time Scat time err flux 
(pe em?) (pc cm?) (s) (s) (s) (Jy) 
312.5 0.01081 j0.010 ja 0.177 
641.8 0.0118 j0.010 6.0 


26 https://www.chime-frb.ca/repeater_catalog 
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Table 5. continued 
687.3 0.00197 j0.0015 1.8 
463.8 0.05407 j0.0047 0.30 
flux_err flux_notes fluence fluence_err fluence_notes sub_num 
(Jy) (Jy ms) (Jy ms) 
0.056 1.48 0.26 0 
2.9 38 16 0 
1.2 1.13 0.69 0 
0.14 - 3.98 0.96 zi 1 
mjd 400 mjd 400 err mjd inf mjd inf err width fitb width fitb err 
(MJD) (MJD) (MJD) (MJD) (s) (s) 
58793.76042828448000 0.00000000000027 . 58793.76032854040000 0.00000000000032 0.00512 0.00045 
59021.048010587610000 0.000000000000053 59021.0478017229 0.0000000024 0.00500 0.00037 
58699.981256104955000 0.000000000000083 58699.981041399835000 0.000000000000088 0.00076 0.00010 
59269.249367994700000 . 0.000000000000096 59269.2492151267 0.0000000018 0.00500 0.00067 
sp_idx sp_idx_err sp_run sp_run_err high freq low_freq 
(MHz) (MHz) 
10.5 2.3 —10.9 2.5 800.2 435.6 
6.9 1.5 —13.4 2.4 800.2 400.2 
14.5 4.0 —26.8 6.8 678.1 400.2 
37 10 —204 54 485.5 400.2 
peak freq chi.sq dof flag frac drift rate R.cc 
(MHz) 
618.5 ed 3.le — 13 
566.6 -1 0.43 
510.8 =! 0.62 
430.7 -1 8.3e — 07 


D. DESCRIPTION OF FIELDS IN CATALOG 


In Table 6, we provide descriptions for each field present in our repeater candidates catalog. An example of these 
entries for four FRBs can be found in Table 5. 


Table 6. Repeater Catalog Field Descriptions 


Column Number Unit Column Name Description 
0 tns name TNS name 
1 previousname Previous name (if applicable) 
2 repeater name Associated repeater name (if applicable) 
3 degrees ral Right ascension (J2000; first likelihood contour; see §3.1) 
4 degrees ra 1 err low Right ascension lower error (see 83.1; first likelihood contour) 
5 degrees ra 1 err up Right ascension upper error (see 83.1; first likelihood contour) 
6 degrees ra 2 Right ascension (J2000; second likelihood contour; see 83.1) 
7 degrees ra 2 err low Right ascension lower error (see 83.1; second likelihood contour) 
8 degrees ra 2 err up Right ascension upper error (see 83.1; second likelihood contour) 
9 ra.notes Notes on right ascension 
10 degrees dec.1 Declination (J2000; first likelihood contour; see 83.1) 
11 degrees dec.1.err low Declination lower error (see 83.1; first likelihood contour) 
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dec. 1 err up 
dec_2 
dec_2_err_low 
dec_2_err_up 
dec_notes 

gl 

gb 

exp_up 
exp_up_err 
exp_up_notes 
exp_low 
exp_low_err 
exp_low_notes 
bonsai.snr 
bonsai.dm 
low ft.90 
up.ft.90 
low ft. 95 
up.ft. 95 
snr fitb 

dm fitb 

dm fitb err 
dm exc. 1.ne2001 


dm. exc.2.ne2001 


dm exc.1 ymwi6 


dm exc.2 ymw16 


bc width 

scat time 
scat time err 
flux 

flux err 

flux notes 
fluence 
fluence err 
fluence notes 


sub num 


mjd.400 


mjd.400.err 


mjd inf 


mjd inf err 


'Table 6. continued 


Declination upper error (see 83.1; first likelihood contour) 
Declination (J2000; second likelihood contour; see 83.1) 
Declination lower error (see 83.1; second likelihood contour) 
Declination upper error (see 83.1; second likelihood contour) 

Notes on declination 

Galactic longitude 

Galactic latitude 

Exposure for upper transit of the source 

Exposure error for upper transit of the source 

Notes on exposure for upper transit of the source 

Exposure for lower transit of the source 

Exposure error for lower transit of the source 

Notes on exposure for lower transit of the source 

Detection S/N 

Detection DM 

Lower limit fluence threshold (9096 confidence) 

Upper limit fluence threshold (9096 confidence) 

Lower limit fluence threshold (9596 confidence) 

Upper limit fluence threshold (9596 confidence) 

S/N determined using the fitting algorithm fitburst 

DM determined using the fitting algorithm fitburst? 

DM error determined using the fitting algorithm fitburst* 

DM excess between DM determined by fitburst and NE2001 as- 
suming the first likelihood sky position of the source 


DM excess between DM determined by fitburst and NE2001 as- 
suming the second likelihood sky position of the source 


DM excess between DM determined by fitburst and YMW16 as- 
suming the first likelihood sky position of the source 


DM excess between DM determined by fitburst and YMW16 as- 
suming the second likelihood sky position of the source 


Boxcar width of the pulse 

Scattering time at 600 MHz? 

Scattering time error? 

Peak flux of the band-average profile (lower limit) 

Flux error 

Notes on the burst flux 

Fluence (lower limit) 

Fluence error 

Notes on the burst fluence 

Sub-burst number (if applicable). If the FRB has only one burst, 


then the sub-burst number is 0. Sub-bursts listed in chronological 
order. 


Time of arrival with reference to 400.1953125 MHz for the specific 
sub-burst. 

Time of arrival error with reference to 400.1953125 MHz for the 
specific sub-burst. 


'Time of arrival with reference to infinite frequency for the specific 
sub-burst. 


'Time of arrival error with reference to infinite frequency for the spe- 
cific sub-burst. 
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S width fitb 
S width fitb err 
Sp.idx 


Sp. idx err 
sp. run 
sis sp run err 
MHz high freq 


MHz low freq 
MHz peak freq 
chi.sq 
dof 
flag frac 


MHz ms ! drift rate 


r.cc 


'Table 6. continued 


Width of sub-burst using fitburst 
Width error of sub-burst using fitburst 
Spectral index for the sub-burst 
Spectral index error for the sub-burst 
Spectral running for the sub-burst 


Spectral running error for the sub-burst 


Highest frequency band of detection for the sub-burst at FWTM 
Lowest frequency band of detection for the sub-burst at FWTM 


Peak frequency of detection for the sub-burst 


Chi-squared from fitburst 

Number of degrees of freedom in fitburst 
Fraction of spectral channels flagged in fitburst 
Drift rate for burst 


Contamination rate (see Section 2.2.2). 


E. 
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Figure 14. Contamination rate (Rec) of all combinations of bursts in clusters of 3+ events. Note that the vertical axes are 
inverted. The gray regions indicate 0.5 < Rec < 5 (silver sample; lighter region) and Rec > 5 (darker region). Potential outlier 
events are marked by red pluses, all other events by black crosses (82.2.2). 
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Figure 15. Continued from Figure 5. 
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Figure 16. Continued from Figure 5. 
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Figure 17. Continued from Figure 5. 
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